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ABSTRACT 

We present a Chandra and HST study of IC 10 X-l, the most luminous X-ray binary in the closest starburst 
galaxy to the Milky Way. Our new hard X-ray observation of X-l confirms that it has an average 0.5-10 keV 
luminosity of 1.5 x 10 38 erg s" 1 , is strongly variable (a factor of «2 in <3 ks), and is spatially coincident 
(within 0."23±0."30) with the Wolf-Rayet (WR) star [MAC92] 17A in IC 10. The spectrum of X-l is best fit 
by a power law with V w 1.8 and a thermal plasma with 4T«1.5 keV, although systematic residuals hint at 
further complexity. Taken together, these facts suggest that X-l may be a black hole belonging to the rare class 
of WR binaries; it is comparable in many ways to Cyg X-3. The Chandra observation also finds evidence for 
extended X-ray emission co-spatial with the large non-thermal radio superbubble surrounding X-l. 
Subject headings: galaxies: individual: IC 10 — X-rays: binaries — stars: Wolf-Rayet — X-rays: ISM 



1. INTRODUCTION 

IC 10 is a metal-poor (Z « 0. 1 5Z^ : Iteaueux et al Jfi979i) . 
barred dwarf irregular in the Local Group. Although ham- 
pered by uncertain reddening corrections, reliable optical and 
infrared distance estimates place IC 10 at 0.6-0.8 Mpc (e .g., 
iSaha et al 17993 iSakai et al 17993 iBorissova et aTll2000t we 
adopt 0.7 Mpc), indicating that it is the nearest starburst 
galaxy to the Milky Way. IC 10 is notabl e for its vigorous star 
formation (^0.03 yr" 1 kpc~ 2 ; e.g.. iHun ter & Gallagher 
ll986HThronson et al.ll990tlHunter et alll993l) and unusually 
large massive (OB) star population, including a galaxy-wide 
surface density of Wolf-Rayet (WR) stars >4 times that ob- 
served in the most active region s of star formation in any other 
Local Group galaxy jMassev & Johnson! fl998t iRover etafl 
1200 It ICrowther et al.ll2003l) . This ongoing star formation in 
IC 10 is expected to result in copious X-ray emission, both 
from compact high-mass X-ray sources and supernova-heated 
gas. The galaxy was observed twice using the ROSAT HRl in 
1996 (70.3 ks combined exposure), although the X-ray flux 
in the HRI band is diminished by a fa ctor of m5 due t o the 
Galactic column of 4.8 x 10 cm" 2 ( St ark et al.l [T992) and 
the e xpected absorption inter nal to IC 10 of «2xl0 21 cm" 2 
(e.g. JYang & SkillmanlH99l . One highly significant source 
was previously detected (hereafter X-l) within the optical ex- 
tent of IC 10 JBrandt et al.lll997l) . This source was found to 
vary by a factor of «3 on ~1 day timescales, with an average 
absorbed HRI flux of 4 x 10" 13 erg cm" 2 s" 1 (or an unabsorbed 
£0.1-2.5 keV ~ 7 x 10 37 erg s" 1 ). The HRI-derived position (5" 
rms) placed X-l in a region with intense star formation as well 
as the most massive H I cloud in IC 10. Notably, X-l was 
found to lie w2" from [MAC92] 17 (V = 21.76), an emission- 
line source identified as a WR star by Ma ssev et al.l (119921) 
and ICrowther et alJ (120031) . and within ss8" of the centroid 
of a non-thermal radio superbubble (w45"/150 pc diameter; 
I Yang & Skillmanl ll993l) . The combination of the high X-ray 
luminosity and strong variability argued that X-l is a powerful 
X-ray binary (containing a neutron star or black hole). 

In this Letter we report on Chandra and HST follow-up 
observations of IC 10 that further characterize the nature of 
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X-l. In particular, these observations more accurately de- 
termine the X-ray luminosity and spectrum of X-l (both of 
which were uncertain since the ROSAT HRI had essentially 
no spectral capability and was affected by the large absorp- 
tion) and confirm the likely association of X-l and the WR 
star [MAC92] 17 (there are several potential optical counter- 
parts to X-l within the 5" radius HRI error circle). In f|2]we 
outline the relevant X-ray and optical observations and anal- 
yses, while in fJ3]we summarize our findings and discuss X-l 
in the broader context of black-hole binaries (BHBs). 

2. OBSERVATIONS AND ANALYSIS 
2.1. Chandra Observations 

IC 10 was observed on 2003 March 12 with ACIS-S (the 
spec troscopic array of the Advanced CCD Imaging Spectrom- 
eter; |Gannii^!uil]20(B|). To limit telemetry saturation in the 
case of high background, only chips S2 and S3 in the ACIS-S 
array were operated, while a 500-pixel subarray was used to 
lower the frame time to 1.7 s and thus limit potential pileup 
of X-l. Events were telemetered in Very Faint (VF) mode, 
and the CCD temperature was -120° C. The optical extent of 
IC 10 (~6.'3x5.'l) is only marginally cropped by the 2-chip, 
500-pixel subarray configuration (16.'8x4.'l). 

Analysis was performed using the CIAO software (v2.3) 
provided by the Chandra X-ray Center (CXC), but also 
with FTOOLS (v5.2) and custom IDL software. The 
data were corrected for the radiation damage sustained by 
the CCDs during the first few months of Chandra oper- 
ations using the Cha rge Transfer In efficiency (CTI) cor- 
rection procedure of Towns lev et alJ (120021) . 2 Following 
CTI-correction, the data were reprocessed with the CIAO 
tool ACIS_PROCESS_EVENTS to remove the standard 0"5 
pixel randomization and to flag potential ACIS background 
events using VF mode screening. We performed standard 
ASCA grade selection in the 0.3-8.0 keV band, excluding all 
bad columns, bad pixels, VF-mode background events, and 
cosmic -ray afterglows. We used only data taken during times 
within the CXC-generated good-time intervals. The back- 
ground varied by <30% during the observation, with the av- 
erage background rate being 2.2 x 10~ 7 counts s" 1 pixel" 1 and 

2 For details see http://www.astro.psu.edu/users/townsley/cti/ 
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7.0 x 10~ 7 counts s _I pixel -1 on chips S2 and S3, respectively 
(in agreement with ACIS quiescent-background calibration 
measurements). 3 The total net exposure time is 29,191 s. 

X-ray sources were iden tified using the CIAO wavelet 
algorithm WAVDETECT JFreeman et all 120021) on the 
0.3-8.0 keV (full), 0.3-2.0 keV (soft), and 2.0-8.0 keV (hard) 
images with a probability threshold of 10~ 6 . In total, 43 
sources were detected. The positions, counts, and spectra 
for these sources were derived using ACIS_EXTRACT (AE; 
iBroos et alJ2002l) with 95% encircled-energy radii [based on 
point spread functions (PSFs) from MKPSF and the PSF li- 
brary]; the adopted X-ray source positions were found using 
the matched-filter method within AE. 

2.2. HST Observations 

To compare the X-ray and optical emission, we observed 
IC 10 with the HST Advanced Camera for Surveys (ACS; 
iFord et al.lfl998h on 2002 October 12. IC 10 was imaged 
along its major axis (ss 50°) using two ACS pointings sep- 
arated by 130". Each pointing consisted of w 0.5-orbit ex- 
posures in the BVI bands (F435W, F555W, F814W). The 
standard pipeline data products were used, although we mo- 
saiced the two pointings ourselves due to a problem with the 
pipeline MULTIDRIZZLE procedure. Optical source posi- 
tions and magnitude s were me asured using SEXTRACTOR 
(v2.2.1: lBCTtin & Arnoutsl 19961) . 
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FIG. 1. — X-ray (dark) and radio (light) contours in the vicinity of X- 
1 overlaid on the HST ACS F814W image. The image shows X-l in re- 
lation to the X-ray and radio emission from the surrounding superbubble. 
Inset: 50-pixel («2"2) close-up of X-l indicating the 0"30 X-ray error cir- 
cle (statistical+systematic) relative to the ACS F814W counterparts. X-l lies 
0"23 from the confirmed WR star rMAC921 1 7A (slightly above and to the 
left of the X-ray centroid; Crowther et al. 2003). 



2.3. Astrometry of Optical and X-ray Images 

Due to the small number of X-ray sources that overlap the 
ACS field of view, we aligned both datasets to the set of 
high-quality ground-based images available from the NOAO 
archive for IC 10 (PI: Massey). Briefly, the NOAO images 
were matched to 55 GSC2.2 sources within 10' of IC 10 
(Aa = 0."20, A<5 = -Of'26), giving an rms scatter of Of'19; this 
alignment was confirmed using 1 1 Tycho 2 stars (0."26 rms). 
Given the optical source density in the IC 10 region and appar- 
ent small, systematic field distortions in the NOAO images, 
we only matched 332 HST sources brighter than B = 23 within 
a 1.'5 region of X-l (Aa = 5."32, AS = 0"27). This gives an 
astro me trie accuracy of 0." 10 rms for the HST positions within 
the vicinity of X-l. Finally, 14 of the 43 significant X-ray 
sources (including X-l) have plausible optical counterparts, 
whereas only two false matches are expected. Aligning the 
X-ray and optical images (Aa = —Of '44, AS = 0."53) yields an 
astrometric accuracy of 0f'28 rms for the X-ray positions. 

2.4. X-ray Analysis 

IC 10 X-l has 4403 counts in the full band and, based 
on our astrometric solutions, it has a J2000 position of 
a =00 h 20 m 29. s 09, 5 = +59°16'51f'95. X-ray contours of X-l 
are shown in Fig. ^ overlaid on the HST ACS F814W im- 
age. The Fig. [0 inset indicates the Chandra positional un- 
certainty of X-l, which is an order of magnitude smaller 
than found previously with the ROSAT HRI. X-l lies 0."23 
from J002029.1+5 91652.1 (an h„ = 2 1.8 WR star identified 
as TMAC921 17 bvlMassev et all 199 2 and [MAC92] 17A by 
iCrowther et al. 120031) . Of'33 from J002029.0+59165 1 .6 {Jab = 
24.4), 0."37 from J002029.1+591651.7 (I AB = 24.0), and 0."42 
from J002029.0+59 165 1 . 8 (an 7^ = 21. 8 OB supergiant iden- 
tified as [MAC92] 17B by Clark et al. 2003). Given the astro- 
metric uncertainty in both the HST and Chandra coordinate 

3 See http://cxc.harvard.edu/contrib/maxim/bg/index.html 
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FIG. 2.— The 0.8-2.0 keV ACIS-S radial profile (small filled circles) of 
X-l, compared to our simple model with and without background (dashed 
and dotted curves, respectively) and the Chandra on-axis PSF (solid curve) 
calculated at 1.5 keV (i.e., closest to where the spectrum of X-l peaks). The 
shaded region indicates the 1 a deviation of the measured profile. Deviations 
from the PSF can be seen out to radii of Ri35 pixels above 0.8-2.0 keV back- 
ground of 0.010 counts pixel -1 . A pixel is equivalent to Of'492. 

frames (i.e., 0."30, much larger than the individual centroid- 
ing errors), X-l is most likely associated with [MAC92] 17A. 
Although we cannot exclude entirely an association with the 
latter sources (or fainter undetected sources) based on the as- 
trometric errors, the probability of X-l lying within Of'23 of a 
WR star by chance is «0.05-0.2%. 4 

In addition to an X-ray point source in Fig. ^ we found 
faint extended X-r ay emission rough l y cent ered on X-l and 
co-spatial with the Yang & Skillman ( 1993) non-thermal ra- 
dio superbubble (see Fig.0. This X-ray emission extends out 
to ^20" in the adaptively smoothed image. Given that the 
X-ray emission appears azimuthally symmetric to ~15", we 
extracted soft and hard-band counts around X-l in 1 -pixel an- 

4 Narrow-band photometric observations of IC 10 have identified 
~ 100 WR stars, of which 26 have been spectroscopically confirmed 
1 Mass ey & ArmandrofH995llRover et all200UlCrowfher et all2003l) . 
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FIG. 3.— X-ray light curve of IC 10 X-1 (dotted) overlaid with the 99.9% 
confidence level Bayesian "blocks" of constant count rate (solid). The events 
are shown in 200 s bins. 



nular bins. Due to the large foreground absorption column of 
IC 10, we used the 0.8-2.0 keV band to reduce background 
where we detect no signal. For comparison, we generated the 
on-axis ACIS-S PSF using CHART, 5 extracting counts in an 
identical manner to the data. The radial surface-brightness 
distributions and la errors for the data are shown in Fig. [2] 
and demonstrate that X-1 lies systematically above the back- 
ground out to radii of ^35 pixels (~ 17") and is clearly ex- 
tended compared to the PSF model. The hard-band profile is 
consistent with a point source and background only. 

We find 175 (129) counts in the full (soft) band between 
radii of 7-35 pixels (note that contamination from X-1 should 
be < 1% at such radii). The fact that the extended emission 
is present almost exclusively in the soft band further supports 
the reality of the extent. An emission model consisting of 
both a point source and a uniform disk with a 17" (60 pc) 
radius convolved with the PSF represents the 0.8-2.0 keV 
surface-brightness distribution well, although there is clearly 
still some residual scatter due to possible clumping. The flux 
ratio of the extended disk component to the point source is 
* 2%. The relative softness of the extended emission is con- 
sistent with a physical origin as supernova-heated hot gas (see 
al so the X-ray spectral analysis below). 

Bran dt et alJ ([1997) found that X-1 varied by a factor of 
«3 during their ROSAT HRI observation but did not have 
the statistics or the time sampling to constrain the variabil- 
ity further. With «4400 counts over a contiguous 29.2 ks, our 
light curve probes the short-term temporal properties better. 
A Kolmogorov-Smirnov test indicates that X-1 is variable at 
the > 99.9% confidence level. To derive the amplitude and 
characteristic timescale of variability, we used the Bayesian 
Block method provided in the contributed CIAO software 
package SITAR. 6 Fig. [5] shows a binned light curve of X-1, 
overplotted with a series of contiguous "blocks" within which 
the count rate is modeled as a constant at the 99.9% confi- 
dence level. X-1 has peak, trough, and average count rates 
of 0.182 ct s~\ 0.088 ct s _1 , and 0.151 ct s _1 , respectively, 
and exhibits rapid variability by as much as a factor of «2 in 
<3000 s. The hardness ratio, however, shows no significant 
spectral variability. Since a significant fraction of X-ray bina- 
ries are X-ray pulsars with pulse periods ranging from 0.01- 
1,000 s (e.g., Whi te et alJfl995T) . we generated power spec- 
tra to search for pul sations using the robust z 2 method (e.g., 
iBuccheri et al.ll9 83). These power spectra also allowed us to 
search for other periodic phenomena (e.g., eclipses or regular 

5 See http://asc.harvard.edu/chart7] 

6 See http://space.mit.edu/CXC/analysis/SITAR/index.html 
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FIG. 4. — X-ray spectrum of the IC 10 X-1 point source, modeled with 
an absorbed power law and thermal plasma with V = 1.83, kT = 1.49 keV, 
and Nh = 6.0 X 10 21 . Note the large residuals around 2.25 keV and 3.65 keV 
suggesting the presence of possible emission lines. 



X-ray bursts). No periodic behavior was found. 

The X-ray spectra were analyzed using XSPEC (vll.2; 
lArnaudl[l99 6). Unless stated otherwise, spectral parameter 
errors are for the 90% confidence level, assuming one param- 
eter of interest. The X-ray fluxes and absorption-corrected 
luminosities for X-1 were calculated from the spectral fitting. 
For the detected X-1 count rates above, pile-up is estimated 
to be ~15%. Thus, spectral analysis of this source must pro- 
ceed with care. To model the spectrum in the presence of such 
pile-up, we used the forwa rd-modeling tool LYNX developed 
at PSU (IChartas et alJ2000l) . The ACIS spectrum of X-1 was 
initially fit with an absorbed power-law model, giving best-fit 
values of N H = (5.1 ±0.4) x 10 21 cm" 2 and T = 1.77 ±0.08 
(x 2 = 97.6 for 82 degrees of freedom). The \ 2 value and 
apparent systematic residuals (bottom panel in Fig. |4ji indi- 
cate the model provides a poor fit to the data. Given that 
we detect a low-level extended component around X-1, we 
added an absorbed thermal plasma component (mekal; e.g., 
iMewe et"afll 19851) to our spectral model. The resulting best- 
fit values are A^h = (6.O+0 Qg) x 10 21 cm" 2 (in agreement with 
the expected Galactic+intrinsic absorption), V = 1 .83^'^, and 
kT = 1.49±8;Ji keV (x 2 = 61.6 for 80 degrees of freedom). 
The abundances were not well constrained and were fixed 
at Z = 0.15Z Q (see While this model provides a sig- 
nificant improvement to both the \ 2 value (>99.99% based 
on the F-test) and the apparent systematic residuals below 
~2 keV, there are still notable residuals around 2—4 keV 
(see middle panel in Fig. @}. These residuals, particularly 
around 2.25 keV and 3.65 keV, are perhaps due to marginally- 
resolved emission lines from an X-ray-photoionized wind 
(e.g., as seen in Cyg X-3; see Sj3}. The 0.5-8.0 keV absorbed 
flux and unabsorbed luminosity of X-1, corrected for pileup, 
are 1.57 x 10~ 12 erg cm" 2 s" 1 and 1.50 x 10 38 erg s" 1 , respec- 
tively. The thermal component contributes «20% of the to- 
tal absorbed flux and w27% of the total unabsorbed luminos- 
ity, which is substantially more than was estimated from fit- 
ting the radial profile. This suggests that either much of this 
component arises from X-1 itself or that the extended com- 
ponent is strongly centrally peaked rather than uniformly dis- 
tributed. We also analyzed the spectrum of the extended X-ray 
component (3."5-17"annulus around X-1). Given the small 
number of counts in this region, w e modeled its spectrum 
using the Cash statistic jCashlll979t) instead of \ 2 - An ab- 
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sorbed mekal model with fixed A/h = 5.9 x 10 21 cm -2 (as found 
for X-l) gives abest-fit temperature of kT = 0.87_iJ'2i keV. 
The 0.5-8.0 keV absorbed flux and unabsorbed luminosity of 
this extended component are 1.73 x 10~ 14 erg cm -2 s" 1 and 
3.24 x 10 36 erg s , respectively. This flux is consistent with 
our radial-profile fitting and supports an origin as hot gas as- 
sociated with the radio superbubble. 

3. DISCUSSION 

This Chandra observation of IC 10 has resulted in improve- 
ments by factors of «15 in terms of spatial resolution and 
astrometric accuracy and wlO in terms of statistics for X-l, 
compared to previous ROSAT HRI exposures, allowing sub- 
stantially better constraints to be placed on the nature of this 
object. The combination of the high luminosity and strong 
variability clearly demonstrates that X-l is a powerful X-ray 
binary, containing a black hole or neutron star. Its likely phys- 
ical association with [MAC92] 17 A implies that the progen- 
itor of X-l must have evolved even more rapidly and been 
more massive than [MAC92] 17 A (WNE stars typically have 
M k, 40-50M(7i and fi ife < 5 Myr; e.g., Maed er & Mevnetl 

1994) . Thus, we expect X-l to be a BHB (altho ugh we note 
that t here are clear counter examples; see, e.g., iKaoer etafl 

1995) . Optical spectroscopic monitoring of X-l to search for 
variability would be one of the best met hods to prove the WR- 
BHB interpretation; Clar k et all ( 120 03 ) find hints of variabil- 
ity in the A4686 feature of [MAC92] 17 A, but longer dedi- 
cated observations are needed to secure this finding. 

The X-ray spectrum of X- 1 is acceptably fit with a T = 1.83 
power law and kT = 1 .49 ke V ther mal plasma (pot entially in- 
dicative a BHB in a hard state; e.g.. lMcClintock & Rem illard 
2003). The residuals of X-l in the 2-4 keV range hint at fur- 
ther spectral complexity, perhaps from emission lines or ab- 
sorption features. Longer Chandra or XMM -Newton observa- 
tions of X- 1 would also be useful to quantify the origin of the 
spectral residuals, as well as to place additional constraints on 
the nature of the X-ray emission (e.g., to rule out further the 
presence of pulsations or bursts). 

The properties of X-l bear several similarities to the 



Galactic X-ray binary Cyg X-3 in terms X-ray luminos- 
ity, spectrum, and vari a bility (e.g., iKitamoto et alJ I1994L 
Liedahl & Paerels 1996; iPredehlet alJ 120001) as well as 
donor type (WNE vs. WN7: Ivan Kerkwiik et alJ I1996L 
Crowther et al. 2003). Moreover, Cyg X-3 has a complex X- 
ray spectrum including numerous emission lines and a 9 keV 
edge that are attributed to an X-ray-photoionized wind; the 
residuals seen from X-l may have a similar origin. Thus 
IC 10 X-l is the only other confirmed example of the short- 
lived WR X-ray binary. 

The discovery of soft extended emission co-spatial with 
the large non-thermal radio superbubble surrounding X-l 
is important for determining the nature of the bubble, al- 
though the poor statistics obtained here do not allow strong 
constraints. Obvious compa risons can be made with 30 
Doradus JDennerl et alJl2001l) . as well as wit h SS433/W50 
JSafi-Harb & Oegelmanl 19971) and IC 342 X-l (goberts et al. 
2003), all of which have similar physical extents and radio/X- 
ray luminosities as the bubble in IC 10. While the radio and 
X-ray emission from the IC 10 bubble is consistent with mul- 
tiple super novae, the superb ubble does lie off the Tj-d rela- 
tion dYang & Skillmanl 19931) . suggesting that something else 
(perhaps X-l) may be powering the expansion. For instance, 
in the case of W50, SS 433 is believed to contribute to the 
expansion. High-resolution radio imaging of IC 10 could be 
used to search for any jets associated with X-l and explore 
whether X- 1 has had any influence on the nonthermal super- 
bubble. 
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